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I. INTRODUCTION 



A system of relativistic particles can be quantized as a constrained system by using the 
auxiliary field formalism (also known as einbein field formalism) to get rid of the square 
root term in the Lagrangian . Applied to the QCD Lagrangian, this technique yields 
Hamiltonians for mesons or baryons H{fii) depending on auxiliary fields /i, 0,|4,0]. The way 
of removing these fields for physical applications is not well defined and can lead to ambigu- 
ities. In practice, these fields, representing for instance the energy density for the particles, 
are finally treated as c-number and determined from the minimum energy conditions 

^=0. (1) 

where E{^i) is an eigenvalue of H{ni). 

When the interaction V between particles does not depend on /ij, the values of yU, which 



minimize the hamiltonian operator (and not E{jji)) are /ij = a/pI + mf. The eigen- 

value equation reduces then to a spinless Salpeter equation with the potential V. So, in this 
case, the auxiliary field formalism can be considered as an approximation of the spinless 
Salpeter equation j^. 

Although handling a relativistic expression for the kinetic energy term is not a serious 
problem with nowadays numerical algorithms, a nonrelativistic expression is always simpler, 
specially for many body calculations. In a sense the auxiliary field formalism simulates a 
relativistic expression with the simplicity of a nonrelativistic one. The price to pay is a 
further minimization on eigenenergies. The values of fii are generally different for various 
radial excited states and this leads to non-orthogonality among them. This is another price 
to pay. An aim of this paper is to study whether these drawbacks are very inconvenient or 
not. 

In Sec. m we consider the case of two identical particles interacting via a potential 
independent of auxiliary fields. We solve the eigenvalue equations for hamiltonians written 
in the auxiliary field formalism and for corresponding hamiltonians in spinless Salpeter 
equations. Then, we compare the eigenvalues and eigenvectors for these two approaches and 
discuss the problem of non-orthogonality. The study is first made for a toy model where 
most of the results are analytical and then in a more realistic one which relies on a precise 
numerical treatment. 
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It is now believed that the color field between a quark and an antiquark in a meson can 
be approximated by a rigid string, carrying both energy and angular momentum. Some 
years ago, a phenomenological model taking into account this particular string dynamics 
has been developed: the relativistic fiux tube model The classical equations of motion 

take the form of a coupled system of nonlinear equations in which a particular quantity is 
introduced, the transverse velocity of the ends of the string, where the particles are located. 
More recently, a lagrangian for a rotating string was derived from the QCD lagrangian 
Adopting a simple straight line ansatz for the string, another coupled system of nonlinear 
equations was obtained, containing auxiliary functions for the particle energy density and 
the string energy density [l^. In this case, the potential does depend on the auxiliary fields. 
In Sec. mil we show that these two approaches are exactly the same, provided all auxiliary 
functions are correctly eliminated from the rotating string hamiltonian. 



II. AUXILIARY FIELDS AND SPINLESS SALPETER EQUATION 

In the formalism of auxiliary fields, the hamiltonian for two particles with the same 
mass m, interacting via a potential V which does not depend on /i, is written j^] = c = 1) 

H^ = ^-^^ + , + V. (2) 
/i 

The philosophy of this formalism is to keep the simplicity of the nonrelativistic form of the 
kinetic energy operator, to get the eigenenergies in terms of auxiliary fields, and to get rid 
of them through a minimization procedure. 

As mentioned, the eigenvalues and eigenvectors of this hamiltonian depends on 

H^\A■nl\^^) = E^l{^^)\A■nl■^J), (3) 

where / is the orbital angular momentum and n (0, 1, . . . ) the vibrational quantum number. 
The physical solutions are obtained by choosing for each eigenstate the value of ^ which 
minimizes its energy. We will note the value of /i minimizing the eigenvalue E^i^n). The 
following notations will be used 

E^i = E^iifXni) < E^iifi) V/i and \A;nl) = \A;nl; f^ni) ■ (4) 
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Applying the Hellmann-Feynman theorem jllj to H , we find 



It follows immediately that 



^lll = {A]nl\p'^ + m'^\A;nl), (6) 



and that 



E^, = 2ii^i+{A;nl\V\A;nl). (7) 

Applied to if the virial theorem yields 2(p^)//i = (r ■ WV). If the potential V is a 
homogeneous function of r of degree a, the formula (jZj) can be simplified. We have then 

= !(l±f!l^„, _ 1^ if F.W = aV. (8) 

a afini 

The scalar product {A;n'l'\A]nl) is proportional to Si^i because of the orthogonality of 
the spherical harmonics. But, for a fixed value of /, two eigenstates with different values of n 
are not orthogonal since they are characterized by different values of /x. So, it is interesting 
to compute the quantity which measures the overlap of two states 

P^,r.i = \{A;n'l\A;nl)\^. (9) 

We have obviously Pnn'i = Pn'ni and Pnni = 1, and we can expect that Pn'ni ~ ^n'n- 

Another way of dealing with auxiliary field is to minimize not the eigenvalues but the 
hamiltonian operator itself and impose the condition dH^ /d^ = 0, fi being considered as 
an operator by its own. In the case of hamiltonian (j21), it is easy to determine that the 
value of the auxiliary field /i which realizes this condition is = a/p ^ + . The minimal 
hamiltonian Q then becomes 



= 2y/p + m2 + V, (10) 

which is simply a spinless Salpeter hamiltonian with the same potential V. We note its 
eigenvalues and eigenstates by 

H'\S;nl) = E^i\S;nl). (11) 
It has been shown in Ref. 0], that we always have E^i > E^^. 
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If the two ways of treating the auxihary field were expected to be very similar, one 
sees that the first method simulates the effect of a relativistic kinetic energy term with 
a nonrelativistic one; the potential is kept unchanged. The differences — give an 
indication on the quality of this procedure. 

Since the minimal energy eigenvalues of are upper bound of the eigenvalues of 
and since these two hamiltonians have in some sense the same physical content, their corre- 
sponding eigenstates must be very similar. By computing the overlap 

Tr,i = \{A;nl\S;nl)\', (12) 

we can expect that Tni ~ 1. 
If we define 

Mni = {S-nlWp + m^\S-nl), (13) 

we can also expect that Mni ~ fini, because of the relation (jH)). In both models, m is some- 
times considered as the current mass, and Mni or fini are then interpreted as the constituent 
mass, which, in this formalism, is state dependent. 

For two massless particles, the relativistic virial theorem [2] yields 2(a/^) = (r • VV). 
If the potential V is a homogeneous function of r of degree a, new relations can be written. 
We have then 

E^, = (l + a)(S;nl\V\S;nl) = '^^^^-^Mni if f-VV = aV and m = 0. (14) 

a 

When the mass m increases, the differences existing between the two approaches and 
tend to vanish and fini ~ rn for all values of quantum numbers. In order to quantify these 
differences, we will study two models containing a confining potential, which is the relevant 
kind of interaction for the hadronic physics. In particular, we will focus our attention on 
the case m = 0, for which we expect the largest discrepancies between the two approaches. 



A. Toy model 

We first consider a "toy model" for which most of the calculations are analytical. Let 
us choose V = kr"^ and m = 0. In this case, is a harmonic oscillator hamiltonian. Its 
solutions are well-known and we have 

E^iifi) = 2 J ^(2n + / + 3/2) + /i. (15) 



The regularized radial part of the wave function is noted 

ut,{r) = {^^kf'^rO^^{{^^kfl'r), (16) 

where 0„;(x) is the usual radial harmonic oscillator wave function and where dimensioned 
factors ensure that the eigenstates \A\ nl; fi) are normalized. 

The values of the parameter fi minimizing the energies are given by 

fini = k'/\2n + l + 3/2f/\ (17) 

and the calculation of the eigenenergies yields 

E^i = 3 k^/\2n + I + 3/2)2/3 = 3 fi^i. (18) 

The relation E:^i = 3 fini can also be obtained directly from formula 

The overlap of two eigenstates can be computed easily. Knowing the values of fini, we 
have 

- I (l^^i^V'° I . (19) 



_2n' + / + 3/2^ 

where the function F„/„;(x) is given, with its properties, in Ref. The quantity Pn'ni is 
independent of k and decreases with increasing values of I and \n — n'\. For instance, we 
have Pi 00 ~ 0.0287, which is the worse case. Considering the states as orthogonal is thus 
always a good approximation. 

In order to solve the corresponding hamiltonian for our toy model, let us consider 
the following nonrelativistic hamiltonian 

:2 



i/^ = |^ + ar. (20) 



Solutions for / = are known [1J|. The eigenenergies are given by 

,2 \ 1/3 



a 



<o = -[^) (21) 

where x„ is the {n + l)th zero of the Airy function Ai(x) {xq ~ —2.338, Xi ~ —4.088, . . . ). 
The regularized and normalized radial part of the eigenvectors is written as 

Ai(x + 



u 



.N 
nO 



(r) = (2z/a)^/6A„((2z/a)^/3^) where A„(x) = ^ . (22) 



Let us note that we have Ai2(g) dq = Ai'^(xr, 
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The hamiltonian is converted into the hamiltonian with V = kr'^ and m = by 



the duahty transformation |p | 



-ar, provided the parameters /c, a, and v are related by 



k 



Due to this relation, we have immediately the property 



E'^, = -{Akf^x^. 



(23) 



In order to compare energies E^q and E^q given respecti vely by relations (fTHjl and (j^ . we 



can use the following approximation for the numbers Xr, 



3vr , , , 
Y(^ + 3/4) 



2/3 



(24) 



With this formula, whose error decreases when n increases, relation ()23j) can be written 

2/3 



i?fo-(|) A:V3(2. + 3/2) 
So we have immediately the ratio (independent of k), 



2/3 



(25) 



£^;^o 3 V 2 



2/3 



0.937. 



(26) 



Numerically, we found that the ratio E^q/ E^q is around 0.944 for n = and tends rapidly 
towards the asymptotical value ()2(j|l as n increases. 

From relation (j^^ and the duality transformation, the regularized function of S-wave 
eigenvectors of can be obtained in the momentum space 



1/3 



P 



(27) 



Consequently, the quantity M„o = {S;nl\^/^\S;nl) can be calculated by the following 
integral = \p\) 

1 



M 



nO 



2V3 



dqqAl{q) = -{4kY/' 
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(28) 



and we find that Mno/fJ^no = E^q/E^q ^ 0.937. The relation E^q = 3 M„o can also be 
obtained directly from formula ()14j) . 

Since the Fourier transform of a harmonic oscillator is also a harmonic oscillator (with a 
phase factor), the overlap T„o can be computed in the momentum space. We obtain 

2i/3 [ ^ - ^ ^ 

^/2n + 3/2 [Jo 



nO 



'^^^^"°\(2n + 3/2)V6 



T7^Mn(2V3,) 



(29) 



which is independent of k. We have checked numerically that Tqo ~ 0.997 and that T„o 
decreases monotonically when n increases. For instance, Tgo ~ 0.523. One sees that, when 
n increases, the eigenvalues for the two formalisms are almost the same (to within 6%), 
whereas the corresponding wave functions can differ appreciably. 



B. Realistic model 



The previous toy model was interesting because most of the quantities can be calculated 
analytically. However, it lacks some physics. First the real confining potential is more alike 
a linear one, and second there exist one gluon exchange contributions which add a coulombic 
term, important at short distance. In order to stick more to the true physical situation, we 
switch, in this section, to a more realistic model. Of course, this study needs a complete 
numerical treatment. 

In the framework of a semirelativistic potential model, a more realistic interaction to 
simulate the dynamics between a quark and an antiquark inside a meson is certainly what 



is called the funnel potential 



V{r) = -'^ + ar with k = 0.5 and a = 0.2 GeVl (30) 

The values chosen here for the parameters k and a can be considered as typical. The 
eigenvalue equations for hamiltonians and with this potential have been numerically 
solved by the Lagrange-mesh method This technique is very accurate and can be 

implemented easily for both nonrelativistic or semirelativistic kinematics. 

In Fig. the quantities AE = {E"^ — E^)/E^ are presented, as a function of the mass 
m, for the three lowest states (IS, 2S, IP) obtained with the potential (jHUj) . As expected, 
AE decreases when the quark mass increases. The corresponding curves for the quantities 
E"^ — E^ present the same profiles. Thus we only present more results for the case of a 
vanishing quark mass which maximizes the error. 

In Table m some energies E^j^ and E:^i are given for n < 3 and / < 2, for the potential (jHOjl 
and m = 0. The difference E:^i — E^i can be quite large, around 100 MeV for the ground 
state. It increases for increasing values of n and decreases for increasing values of /. The 
quantity M„; is also given in this Table; it increases with n and /. The parameter ^ni can 
differ from Mni by several tens of MeV, but the curves E^i^fi) always present a fiat minimum 
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around Uni- So it is not necessary to obtain a precise value of Uni to obtain a good value of 

The overlap Pnn'i is given in Table ITTl for n,n' < 3 and for / < 1, for the potential 
and m = 0. Except in the case \n — n'\ = 1, the overlap is very small. So the eigenstates 
\A] nl) can be considered as quasi-orthogonal. This conclusion is the same than in the toy 
model, and we can hope that it is a universal conclusion, valid whatever the potential V . 

In table UTTl the overlap Tni is presented for n < 3 and / < 2, for the potential (jHUI 
and m = 0. It increases for increasing values of / and decreases for increasing values of n, 
and present similar features to the toy model. For small values of the vibrational quantum 
number, the eigenvectors \A\nl) and \S\nl) are rather similar. 

III. ROTATING QCD STRING 

Up to now, the confining term was put by hand. There exist more serious explanations of 
this term, based on QCD arguments. In particular, one can consider a string with an energy 
density between the quark and the antiquark. This string is itself a dynamical object. Here 
we examine two such approaches. 

We first quickly present the relativistic flux tube model, which is essentially phenomeno- 
logical, and then the rotating string model within the formalism of auxiliary fields, which 
has a more firm foundation on QCD theory. Lastly, we show that these two models are 
completely equivalent. 

A. Phenomenological model 

In the simplest version of the relativistic flux tube model Jfj, a quark and an antiquark, 
with the same mass m, move being attached with a rigid flux tube, assumed to be linear 
with a uniform constant energy density a. The system rotates in a plane with a constant 
angular velocity around the center of mass, which is assumed to be stationary. Denoting 



total angular momentum, v± the transverse velocity of the quark or the antiquark (relative 
to the string direction), 7± = (1 — and U{r) a potential used to describe dynamical 

effects coming from mechanisms other than the flux tube, the classical equations of motion 



r the distance between the two particles, pr the radial momentum 
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of the system are given by 

L 



where 



^ v^-f^^/j^f+m? + ar f{vA_), (31) 
if = 27^V^^+^+ar^^^^^ + [/(r), (32) 



f{x) = -^-^ ^arcsinx — xVl — x^j , (33) 



is a very important function in the formahsm. These equations have been generahzed in 



the case of asymmetrical systems 



|. Numerical solutions of the quantized versions of these 



equations have been obtained 0, S Q, 13 (^^^ ^ practical calculation L must be replaced 
by + 1)). Equations (jSH)-(IS21) are coupled nonlinear equations. In practice, the value 
vj_{m, a, L; pr, r) is extracted from Eq. (j!?T|) and injected into Eq. (jS2I), giving a L dependent 
hamiltonian H{m,a, L]pr,r), which is diagonalized afterwards. 

When the relativistic flux tube hamiltonian is supplemented by appropriate potentials 
(Coulomb-like, instanton induced effects) and when it is assumed that each extremity of 
the flux tube can give a constant energy contribution, rather good meson spectra can be 
obtained hAlm. 



B. Auxiliary fields for a rotating string 

Starting from the QCD lagrangian, a lagrangian for a meson can be derived taking into 
account the dynamical degrees of freedom of the string For a system containing a quark 
and an antiquark with the same mass, with the hypothesis of a straight line configuration for 
the minimal string, and introducing an auxiliary field fi for the energy density of the quarks 
and another auxiliary field u for the energy density of the string, the following classical 
hamiltonian can be obtained 



Mr) + 2 J [P-^j i^{P,r)dp 

where U{r) is a potential describing dynamical effects coming from mechanisms other than 
the string. In this formula, r is the common proper time of the two particles and /3 is a 
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coordinate along the string. Within this formahsm, the two auxihary functions /i and u are 
to be varied and to be found from the minimum of H. 

Let us first remark that, if we neglect the integral in the denominator of the term de- 
pending on LF' or if we put L = 0, then H is minimum for v = ar, independent of (3. In this 
case, because of the relation H — - = p^, H reduces simply to with V{r) = ar + U{r). 
The linear confining term appears naturally in the formalism. A further minimization on fi, 
as demonstrated in second section, gives rise to with the same potential V{r). 

This result is interesting and simple, but indeed one can make an exact minimization on 
the field u. A straightforwards calculation shows that the minimization condition dH/ du = 
is fulfilled if the field u is set to uq with 

uoiP) = =, (35) 



where y is to be found from the transcendental equation 

A = /(!/) + -• (36) 
ar^ ar 

The function / defined previously by relation ()33|1 appears again curiously (for a practical 



calculation L must also be replaced by + 1)). Using expression ()35|1 . one obtains from 
relation ^ 

pl + m^ ,^ 2n arcsiny , 
H = ^ + fx{l + y^) + ar - + U{r). (37) 

/i y I 

WKB solutions for the system (jHBj) and ()37p are obtained in Ref. [10| for m = and 
U (r) = 0. Regge trajectories are computed in agreement with experimental data. 



C. Equivalence between the two string models 

In principle the way of solving the system ()36|) - ()37|) is the following. One extracts 
y(a, L;r;/i) from Eq. p6|) and inject it into Eq. (j37jl leading to H{m,a, L;pr,r; fi). This 
hamiltonian is then diagonalized and, for each state, the optimal value of /i must be deter- 
mined to get the physical eigenenergies. In view of the form for various expressions, this 
procedure seems hopeless. Fortunately this is not so and one can go one step further and 
perform the minimization on /i. 
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The partial derivative of Eq. (jHBj) with respect to /x gives the following relation 

0, (38) 



HU 1 arcsin y 



in which the expression dy/dfj, appears. Now imposing dH/dfj, = 0, we obtain from Eq. (|n7|) 



+ ^ 2 9y 

+ 1 + 2/ + 



a r arcsin 



2 



(39) 



Combining the two last formulas, we obtain 



Using this relation, it can be verified that Eqs. ()H8|1 and are well defined for y G [0, 1]. 
Equation (jinjl giving fi{y) and the presence of the function / (jHHj) in Eq. are the clues 



for the identification of both methods. Replacing y by v± and /i by 7_la/ p.^ + transforms 
Eq. (jnni) into Eq. and Eq. (jHTj) into Eq. proving thus the complete equivalence 

between both approaches. In passing we have now a clear idea of the physical content for 
the auxiliary field fi and of the mysterious variable y appearing in the theory. 

IV. CONCLUSION 

Effective hamiltonians for hadron dynamics can be derived from the original QCD la- 



grangian. These hami, 
quark energy density 



tonians can depend on auxiliary fields, representing for instance the 
^0]. When the interaction does not depend on these fields (exactly 
or by approximation), it is possible to reduce the effective hamiltonian to a semirelativistic 
one (spinless Salpeter). 

In the first part of this paper, we have compared the eigenvalues and eigenvectors for 
hamiltonians written in the auxiliary field formalism and for the corresponding spinless 
Salpeter hamiltonians . We have only considered the case of two particles with the same 
mass and two different confining interaction: A "toy" quadratic potential and a more realistic 
funnel potential. We have shown that the eigenvalues of are close to the eigenvalues of 
, but the relative differences can exceed 10%. The overlap of corresponding eigenvectors 
for the two approaches is generally close to unity for ground states, but it deteriorates 
when the vibrational quantum number increases. It appears that, if precise calculations 
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are desired, it is better to solve directly a spinless Salpeter equations: Accurate numerical 



techniques exist for two-body [,18i], and many-body 



22| problems. Hamiltonians with 



auxiliary fields can be useful if one searches to obtain analytical results, as it is the case 



radial excited states being quasi-orthogonal. 

The phenomenological semirelativistic fiux tube hamiltonian has been developed in order 
to take into account the string dynamics in a meson 0, 1^. Another model derived from 
the QCD lagrangian, the rotating string, has been more recently developed [l^. But the 
rotating string hamiltonian depends on auxiliary fields. 

In the second part of this paper, we have shown that the classical equations of motion of 
the rotating string reduce exactly to the classical equations of motion of the semirelativistic 
fiux tube, provided all auxiliary fields are correctly eliminated. The complete equivalence 
of both models has two virtues: to provide a natural interpretation of the auxiliary field in 
term of the transverse velocity and to reinforce the relevance of the relativistic fiux tube 
model (which was considered up to now as a phenomenological one), since it is equivalent 
to a model based on firm QCD grounds. 

The set of coupled equations for these two models are very complicated to solve numer- 
ically . The difficulties are essentially the same for both approaches, with the 
supplementary necessity to perform a minimization within the rotating string formalism in 
order to determine the value of the auxiliary field for the quark energy density. It is then 
more interesting to work directly with the semirelativistic fiux tube model. 
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TABLE I: Values in GeV for the quantities E^^, E^^ and Mni (see Sec. as a function of the 



quantum numbers n and /, for the potential (|.SOj) and m = 0. 

n 1 2 3 

/ = E^i 1.197 1.912 2.461 2.914 

E^i 1.294 2.109 2.700 3.185 

Mr,i 0.422 0.588 0.712 0.817 

1 = 1 E^i 1.759 2.316 2.782 3.187 

E^i 1.826 2.472 2.990 3.434 

Mni 0.508 0.644 0.757 0.855 

1 = 2 E^i 2.170 2.643 3.056 3.426 

E^i 2.225 2.780 3.248 3.659 

Mni 0.594 0.711 0.813 0.904 



TABLE II: Values of Pnn'i (see Sec.^, as a function of the quantum numbers n, n', and /, for the 
potential ()30|) and m = 0. Values for Z = (/ = 1) are given in the upper-right (lower-left) triangle 
of the Table {Pnn'i = Pn'ni and P„„i = 1). 



Pn n' I 


n = 


1 


2 3 


n' = 


1 


0.022 


1.6 10"^ 2.4 10"^ 


1 


0.017 


1 


0.022 5.7 


2 


2.0 lO-'^ 


0.023 


1 0.028 


3 


5.9 10-5 


2.5 lO-'^ 


0.019 1 



TABLE III: Values of T^i (see Sec. as a function of the quantum numbers n and /, for the 
potential H30() and m = 0. 



Tnl 


n=0 12 3 


I = 


0.9837 0.9639 0.9389 0.8982 


1 


0.9908 0.9722 0.9406 0.9005 


2 


0.9925 0.9746 0.9443 0.9008 
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FIG. 1: AE = {E^ - E^)/E^ (see SecHH as a function of the mass m, for IS, 2S, and IP states 
with potential (j3fl|) . 
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